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FOREWORD 


This Technical Report was prepared by RCA Laboratories, Princeton, 

New Jersey, under Contract No. F30602-72-C-0486. It describes work 
performed from April 1973 to August 15 1973 in the Communications Research 
Laboratory, Dr. K. H. Powers, Director. The principal investigator and 

project scientist was Dr. D. A. de Wolf. 

The report was submitted by the author in August 1973. Submission 
of the report does not constitute Air Force approval of the report s 
findings or conclusions. It is submitted only for the exchange and 
stimulation of ideas. 

The Air Fore'' Program Monitor is Capt. Darryl P. Greenwood. 


SUMMARY AND ABSTRACT 


This Final Report of Contract No. F30602-72-C-0486 regarding the 
behavior of laser beams in turbulent air extends results summarized in 
a previous report (RADC-TR-73-162). It is shown that amplitude scintil¬ 
lation is unimportant in focused beams, and dominant in collimated and 
diverging beams. The limitations of an existing numerical procedure for 
computing these are investigated. Intermittency of turbulence is rede¬ 
fine! in terms of large-scale variations in the refractive-index struc¬ 
ture constant and the dependei.ee upon averaging time is determined. It 
is concluded that intermittency is essentially ar. extra randomness dur¬ 
ing operation of the laser for short periods (of several minutes or 

Significant fluctuation rates are defined for angle of arrival, 
and amplitude. They differ greatly. 
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1. INTRODUCTION 


This is the Final Reprrt of the RCA effort on the effect of turbu¬ 
lence instabilities on laser propagation under contract no. 

F30602-72-C-0486 to Rome Air Development Center. Previous reports under 
this and under the preceding effort (contract No. F30602-71-C-0356) are 
summarized [l] in the list of references, and they will be referred to 
as TRI to TRVIII in chronological order. 

Aside from previously reported developments in TRI-TRVII, also sum¬ 
marized in TRVTII, this work has led to the following developments: In 
section 3, a generalization of the Pythagorean sum of focal-spot area 
and induced broadening by turbulent air is given for any image plane of 
a gaussian beam (diverging, collimated, or focusing). We prove in sec¬ 
tion 4 that amplitude scintillation in the focal plane of a focused beam 
is negligible, and therefore restrict the analytical development of 
laser-beam irradiance to the Rytov approximation in that case. Except 
for exceptional situations, the preceding saturation-regime calculations 
for plane and spherical waves suffice for non-focused beams. Herrmann 
and Bradley's numerical procedure is investigated in section 5 and it 
is shown that this procedure cannot predict large amplitude fluctuations. 
The problem of intermittency is discussed in section 6. The power spec¬ 
trum of at lrw wave numbers is its determining factor, and we con¬ 
clude that intermittency must be regarded as an unpredictable random 
effect governed by long-term statistics in view of the fact that the above- 
mentioned power spectrum is highly specific to terrain and local weather 
conditions. The idea that the angle-of-arrival spectrum can yield micro- 
sca'e information unfortunately must be abandoned (section 7). However, 
in se'tion 8 we show that important fluctuation rates of angle of arrival, 
amplitude, and phase can be defined and computed. Section 2 is an intro¬ 
duction to the problems defined in the whole effort, and to the method 
of solution: it serves as an aid to the reader who does not wish to im¬ 
merse himself in all the details of previous reports and work. 
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2. INTRODUCTION TO THE LASER-BEAM PROPAGATION PROBLEM 

The work of the contract effort of which this is the Final Report 
is centered around the problem of forming an image in the plane z = L 
of a laser beam originating at z = 0. The beam must travel through 
(possibly turbulent) air. One might distinguish between three aspects 
of the problem: 

(i) The optics of imaging a source at z = 0 in the plane z = L in the 
idealized case of free space between source and image. Thus, one 
might consider what optical functions (e.g., mutual transfer func¬ 
tion, point spread function, etc.) give optimal information about 
the source or the imaged object. 

(ii) The effect upon the imaging process of the air through which the 
beam travels from source to image plane. Air can absorb energy 
from laser beams; aerosols and other particulate matter can scat¬ 
ter energy out of the beam, and inhomogeneities in the gaseous 
properties (e.g., in temperature, density, water-vapor content) 
can give rise to deflection of rays. All of these effects influ¬ 
ence the shape of the image. 

(iii) The nature of the source or of the object to be imaged. In our 

work, the source is nearly always the output of a laser, and hence 
it is the shape of the transmitted laser beam that plays the role 
of a separate input to the problem. 

The distinction has been overemphasized. Figure 1 illustrates the 
most important aspect of the problem: the interaction of light rays with 
air. A host of effects occuis even when the intensity is not extremely 
high (we have not sketched in an equally large host of other effects such 
as air breakdown which occur at extremely high intensities). It is im¬ 
portant to realize, however, that the beam quality mentioned in (iii), 
i.e., whether one uses focused beams, or whether the beam is gaussian 
with or without truncation, has a definite effect upon the beam quality. 
For example, a highly focused beam will create higher intensities at the 
focal spot, and one might well expect intensity-dependent distortions 
such as thermal bending and blooming of these to occur more readily than 
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of diverging beams. We shall see that the nature of turbulence-induced 
distortions is influenced importantly by beam shape. Also, the distor¬ 
tion of the image will or will not be important depending upon what as¬ 
pects of the image are important in applications (tracking, weaponry, 
illumination, etc.), fin d consequently the study of which optical parameter 
should be considered [as indicated under (i)] is also of some importance. 

A unified mathematical description of all of these aspects is based 
upon the concept of the electric field E(r) at a location t = (p,L) in 
the image plane. Under a number of not overly restrictive conditions 
(to be discussed later), it is given by 


S(r) = L _1 exp[ik(L + p 2 /2L)1 

x ik C. 2 a 

2^ J d P 1 B (r,r 1 )U o (r 1 )exp[ 


ikPi 


2* ' ( 1 ■ i' -ikp’PjA 1 , 


( 2 . 1 ) 


where r^ - (ppO) is a point on the aperture plane. The other symbols 
are: 

k = bi/c = 2 it/ A : wavenumber of the laser radiation (u) = frequency, 

A = wavelength). 

L : propagation distance 

R : radius of curvature of the phase front at z = 0 ( R > 0 

for a focused beam). 

u 0 ( r l) : aperture or pupil function 

B ( r J r ^) : the ratio of the electric field at r due co a point 

^ * 
source at r in air to that in free space. 

The optical imaging properties we desire depend upon the use we make 

of Eq. (2.1). We may be interested in its phase properties (e.g., for 

angular deviations) or its intensity pattern (e.g., for illumination or 

for weapon uses). The nature of the beam is contained in the choice of 

R and of U (r.). For an untruncated gaussian beam, one has 
^ 2 2 2 

V r l> = ex P("P^ / r c )> where r Q is an effective radius-variance pa¬ 
rameter. For some applications it is more appropriate to take u (? ) = 

•k o 1 

const, for < r^, and U^(r^) = 0 for > r^: thjs choice yields the 

well-known Airy diffraction pattern in the focal plane (R = L) of a 
focused beam. The most crucial aspect of the program - the interaction 
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of laser beams with (turbulent) air - is hidden in the function B(r,r^): 
the response to a spherically symmetric point source in turbulent air. 

The distinction (i)-(in) has been useful to the extent that it 
allows us to separate beam and imaging problems from the air-wave inter¬ 
action. It is a great simplification to study this interaction for only 
a spherical-wave source. Nevertheless, we will not adhere rigidly to this 
way of separating the problems; for some cases it will be easier to look 
at the interaction of air with beam waves. 

Consider the basic interaction problem. A quasi-monochromatic point 

source at r, =0 yields an electric field E that satisfies the equation, 

1 1 o 

V 2 E + k 2 E = -t(r) , (2.2) 

o o 

where b(r) is a Dirac-delta function. The free-space solution for E 

+ -l 

is (aside from some constant factors involving x) E (r) ^ r exp(ikr). 

In air we have 

V 2 E + k 2 £ (r, t)E = -6(r) , (2.3) 

where e(r,t) is the dielectric permittivity. The equation is approxi¬ 
mate to the extent that it ignores terms of order de ( r,t)/u>dt as well 
as terms of order k | Ve(r, t) |/e (r, t). Optical spatial and temporal fre¬ 
quencies are so high that such approximations are warranted. Workers 
in optics often prefer to use the refractive index n(r,t) by substituting 
C(r,t) = n 2 (r,t) in Eq. (2.3). Obviously, the source function - a delta 
function in this case - can be eliminated by subtraction of Eq. (2.2) 
from Eq. (2.3), and the electric field E can be expressed in terms of 

E . If the source is at r , then both E and E are functions of r and 
*° ° 0 * 

r . The previously mentioned normalized field B(r,r Q ) is simply the ratio 

E(r,r )/E (r,r ). 
o o o 

It is possible to simplify Eqs. (2.2) and (2.3) by utilizing the ex¬ 
tremely small size of optical wavelengths With respect to scale sizes of 
variations in e(r,t). Under the so-called paraxial approximations, a 
parabolic equation for B(r,t) can be developed, namely 

V 2 t B + 2ik(r* V - ^)B + k 2 6fB = 0 (2.4) 
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Jk 

where ^ is the gradient-operator component transverse to the radial di¬ 
rection, and 6e is the deviation of e(r,t) from unity. This equation, 
and its plane-wave corollary, which differs from Eq. (2.4) to the extent 
that (r«V - r ) is replaced bv b/'Az, have been studied extensively by 
workers in propagation of waves through air. 

The permittivity deviation bc(r,t) in Eq. (2.4) is determined by 
the refractive index n(r,t). In air, at optical frequencies, the latter 
is given by 


n (r, t) = 1 + 


10‘ 6 (1 


-2 


+ 0.0075A ) x 77.6 p(r,t)/T(r,t) 


(2.5) 


where pressure p is given in millibars, temperature T in degrees Kelvin, 
and wavelength A in urn. The main cause of variations in n(r,t) are those 
in the temperature. The so-called non-linear variations in T due to air 
heating by the laser beam itself are not the topic of this contract effort. 
Rather, the irregular fluctuations in T caused by turbulence - particularly 
of the boundary-layer type - are the subject of this work. Thus an at¬ 
mospheric model for turbulence, possibly homogeneous turbulence, at 
altitude x is required. Tatarski [2] has elaborated on the connection 
between temperature and refractive-index statistics. Aside from constant 
factors one is given by the other. Consider therefore the Fourier trans¬ 
form of the covariance of the dielectric permittivity 


<KK) = 



( 2 . 6 ) 


* v -G r is a point at altitude <. The atmospheric-turbulence model spec¬ 
ifies a spectrum, 

<KK) = 32n 3 x 0.033 C 2 (x)(K 2 +L " 2 ) -U / 6 exp(-K 2 /< 2 ) , (2.7) 

n o m 


where macroscale L is presumably a function of x, microscale is re- 

° - 12 ° 
lated to K by the relationship ( = 5.92k , and C (x) is the refrac- 

m o m * n 

tive-index structure constant related to its sea-level value [3] by the 
relationship 


C n 2(x) = C n 2 (x o )(T o/ T) 2 (X o / x)2/3(1 + X /V ’ 2/3 
x exp[-2(x-x Q )/h] , 


( 2 . 8 ) 
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where x - 1 m and T are sea-level values of effective altitude and tempera- 
o o 

ture, and • 71 is the Monin-Obukhov length (H ■+ °° at dawn and dusk, and 
s s 

?. 'v 1,5 m at midday for well-developed turbulence). The structure func- 
s 2 

tion C (x) can be obtained from temperatrre-variation measurements, and 
n -1 

the microscale k , although harder to estimate, enters into the calcu- 
m ’ 

lations in such a manner that an error of a factor two or so in its mag¬ 
nitude is not serious. The macroscale L varies more but it often drops 

o 

out of the end result. 

These considerations, while not complete, give the essenti'il inputs 
to the calculations made undej the auspices of this program. A previous 
report TRVIII [l] gives a summary of results obtained up to February 1973. 
Ensuing work will be discussed from here on. 
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3. GENERALIZATION OF THE FOCAL-SPOT BROADENING RESULT 


Laser beams cannot focus into a point image because diffraction at 

the edg°s of the aperture through an angle ~ l/kr Q yields a blurring 

at distance L of order L0j ~ L/kr (remember that r is a measure of the 

d ' o o 

aperture radius). In turbulent air, the angle 0 is randomly modified 

d 2 1/3 

into 0, f 60, and 60 has a variance of the order of C L K ' in homo- 
d ’ n m 

geneously turbulent air (e.g., for horizontal propagation). The focal 
spot is broadened in area on the average. 

In TRI, we defined an aperture radius for arbitrary aperture func¬ 
tion U (r. ) - see Eq. (2.1) - through tie definition 
o 1 

r ~ 2 = - fd 2 P U V 2 u / /VpU U (3.1) 

o J 1 o T oj ^ o o 


This definition corresponds to the one-sigma width of an untruncated gaussian 
aperture function, U Q (r^) = exp(-r^ / r Q )• The integration and the gra¬ 
dient operator in Eq. (3.1) pertain to the two variables in the z = 0 
plane. An irradiance-weighted radius r^ was defined in the plane z = L 
by means of an integral, 

r L 2 = /d 2 0 p 2 <K?,D>//d 2 p< I(P.L) > (3 ’ 2) 

and a free-space equivalent, r , is defined similarly by utilizing 
I Q (p,L) instead of I(p,L) . Actually, some first-moment integrals 
need to be subtracted, but they are zero for axially symmetric U Q (p^) 
and for homogeneous turbulence. With these definitions, we obtained a 
very general result, 



2 2 
r T L + 12n x 
Lo 


0.033L 2 < l ^ 3 J L 


ds 


m 


2 2 

s u (s) 
n 


(3.3) 


.for an axial path from s = 0 to s = L. For horizontal propagation, 

r 2 -r 2 = r 2 ~ 1.3C 2 L 3 < 1 / 3 . Conditions for validity were given in 
L o “ LB n m 

TRI, and comparisons with other foc-:l-area definitions were made in TRIV. 

The above definition requires care because r^ Q , hence also r^, become 

infinite for sharply truncated beams. This difficulty is easily circum- 

2 2-2 

vented in practice by processing for r^ -r Lo = r^ (which is always 
finite) but the definition weights the side-lobe structure unduly. On 
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Che other hand, r^g^ * ias ^ unc ^ ament ' a ^ significance as the "displaced area 
of the focal spot due to average turbulence broadening. It is given in 
Eq. (3.3) for the focal pianc. We will now generalize it to any image 

plane. 

Schmeltzer [4] has given the following form for the electric field 
of a gaussian laser mode in the plane z and close to the axis: 

2 

r.1 . • i. l r\ / _ _• 1 _ ^ \ 1 

} 


E (r) = — exp[ikz + ikp 2 /2(z-ikw 2 ) 

° z-ikw 


(3.4) 


-2 -2 -1 
w = r + ikR 
o 

This expression is subject to 


! 


: radiation-zone . m'ition 
: near field of t aperture 

: sagittal conditions 


(i) kz » 1 

(ii) z « kr 2 

3 ° 4 

(iii) z >> kr 

3 ° 

(iv) z >> kp^ 

Conditions (i)-(iv) are sufficiently liberal to include all of the image 
in most cases when the propagation distance is at least a few hundred 
yards. It can be seen from the absolute square of E Q (r), after straight¬ 
forward algebraic manipulation, that 

I (?) = (kr 2 /z) 2 [l + (l-z/R) 2 (kr Q 2 /z) 2 ] 1 
O O 

x exp{-(kr Q p/z) 2 /[l + (1-z/R) (kr Q / z ) H 

2 2 

Thus, the free-space intensity at r = (p,z) is s^ill a gaussian exp(- /r Q ) 
where r is another constant radius determined from Eq. f3.5). The deri¬ 
vation of Eq. (3.3) from Eq. (3.2) is practically the rame, except that ^ 
we obtain ^(z) instead of r^ 2 to get the result; r L ~ = r Lo (z) + r^ B 
with 


r Lo 2(z) = r 0 2 ^ 1 * i / R)2 + ( z / kr 0 2)2] 


(3.6) 


where r 2 is given by the second term of Eq. (3.3). This result is 

IwB *i 

valid at any plane z = const [subject to the conditions (i)-(iv)J. In 
that sense, it is more general than Eq. (3.3), but it is less general in 
the sense that a gaussian mode has been assumed for the aperture funcli.cn 
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U (r). Note that Eq. (3.4) results from Eq. (2.1) upon inserting 
o ^ 2 2 

U (r^ = exp(-p x /r Q ) into the latter equation. 

At z = L = R, we obtain from Eq. (3.6) the previous result. For 

collimated beams (R -*• °°), we note that diffract .on is usually a small 

edge effect because r T 2 (z) = r 2 + (z/kr ) 2 in that case, and because 
° Lo o o 

z/kr 2 << 1 so that r, (z) = r . Our previous studies (see Fig. 3 of 
' o Lo o 

TRVIII) have shown tlat r << r for most cases. Therefore, in the case 
of collimated beams - and certainly for diverging beams - the beam¬ 
broadening effect is less important. In those cases, amplitude scintil¬ 
lations are the dominant intensity effect. This immediately poses the 
question: do amplitude scintillations then also plav an important role 
in the focal plane of a focused beam? Evidently, from Eqs. (3.3) and 
the preceding qualitative explanation on the basis of refractive bending 
of rays superimposed upcu the tree-space diffraction effect, phase effects 
are important. This question is answered in the next section. 
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4. BEAM WAVE IRRADIANCE AND PHASE IN TURBULENT AIR 


We shall deal here with amplitude scirtillation in the focal plane. 
Previously, we developed a comprehensive Lheory for plane waves (TRV and 
TRVI) and frr spherical waves (TRVIa). These theories are important for 
collimated and diverging beams respectively wher (see previous sectior) 
we have shown that refractive (i.e., phase) effects are only important 
at the edge of the beam. W? shall derive a Rylov approximation for 
gaussian-beam waves analogous to the development in sections 1 and 2 of 
TRV for plane waves and in s ctior 2 of TRVII for spherical waves. It 
appears useful to give a very abbreviated review of the steps of the 
derivation. First of all we note that the normalized field B can be written 
as a sum, 


B = 1 + B x + B 2 + B 3 + - 


B „ " /i 3 r l G(t,# 1 )«erf l )--yj 3 r ii G(f ii . 1 ,t ii )«e:(r i ) , 


(4.1) 


where G(r,r 3 ) is a Green's function for beam waves (all trivial factors 
can be absorbed in G). The next step consists of utilizing the small- 

^ mM 

angle scattering properties of fie which imply that r m ^" r m makes only 
a small angle with the z axis. By means of a stationary-phase analysis 
(or a steepest-descent method one can perform all the transverse inte¬ 
grations in (4.1) to obtain formally something like 

z 


/ n-1 

dz 1 f(z 1 )6e(0,z 1 )- f dz f(« )«e ( P,« )g( V ~« ) , 


(4.2) 


where f(z ) is a- function resulting from (4.1) after the transverse in- 
m 

tegrations. The problem is that there is also an inseparable exponential 

function of z.-z which we denote by g(z,-z ) that prevents us from 

In In 

separating the integrations over dz^ (1 £m £ n) in general. The Rytov 

approximation is obtained bv finding the conditions under which 

g(z-z ) = 1. In that case (4.2) reduces to 

1 n 


rl 

B n - ~T [J dzf(z) fie(0,z)] n 


(.4.3) 
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T 



(4.4) 


and upon summation over n as in (4.1) we obtain 
L 

B = exp[ J' dzf (z)fe(0,z)] 
o 

We have abbreviated the actual procedure too much because it really has 
been carried out for moments of B rather than for B - thus yielding energy 
conservation which Eq. (4.4) does not yield - and we have omitted some 
intermediate extra steps such as utilizing transverse Fourier transforms 
and many statistical properties of $£, but the many details of the analysis 
are given in the above references (and exhaustively in TRIII). 

Let us now give the details of this procedure for gaussian laser 
modes given by Eq. (3.4). Referring to section 2 of TRIII, we note that 
the beam-wave Green's function G^r,^) is given by 

= G p ( -l )E o ( V/ E o^) (4.5) 


where E is given by Eq. (3.4), and G (r-r, ) by Eq. (2.7) of TRIII. This 

u pi 

Green's function is now gernune in Eq. (3.1) of TRIII. In order to per¬ 
form the d integrations ir that equation, we perform a stationary- 
phase (steepest-descenf ) analysis on the integrand as in section 6 of 
TRIII. The stationary-phast function in the exponent of the integrand 
Is 


V - k[AP m ! + Az V / 2 + p 2 /2Z - p . 2 /2Z J-K.p 
n m m m ' m M m-1 ' m-1 mm 


Z - [ z - i kw “ ] 
m m 


(4.6) 


-2 . 
w - 



+ ikR 


The points of stationary phase are found by setting the gradient of 

^ m 

with respect to p^ equal to zero. A result very close t •> the spherical- 

wave result of TRVII [Eq. (2)] i found; 


A V Ar m = -%/ k + (4-7) 

The notation is consistent with previous work, and the interpreta¬ 
tion of Eq. (4.7) is practically identical to that of Eq. (2) of TRVII. 
Again, small-angle scattering determines the physics, and small-angle 


1? 



approximations in Eq. (4.6) are permissible, yielding after some rearrange* 
ment • 


Y 

m 


. Z Z . /p , 
k m m-1 / m-1 

22 ,-Z Iz . 

m-1 m \ m-1 




(4.8) 


This result is al 


From it, we derive, as 

n ik h 

B - .n, -2 J 


J-l 


8 t 


so quite parallel to that for spherical waves, 
in section 2 of ’’RVII but now for beam waves: 



x n expt-iK m 2 (L-, ai )2 ! /2kZ o ] e xp[-iK ni .q iitfl 'a- Zm )/k] 


m : 


=1 


Q - ' = V (Z./Z )K 

Vi Mi J ° J 


(4.9) 


2 

Z = L-ikw 
o 


We obtain the normalized field B from (4.9) by summing over all in¬ 
teger n from n = 0 on. Here too, as in section 17 of TRIIT- the Rytov 
approximation (modified automatically so that energy is conserved) fol¬ 
lows if the second exponential in (4.9) is very close to unity. In that 
case, the parameter can be shown to be [in analogy to Eq. (17.10) 

of TRIII without the last exponential factor]: 


<6x 2 > 



dKK < J'(K)e" K n ( z )( L-z )/ k |l-cos[K 2 £(z)(L-z)/k]j 

(4.10) 


where we have utilized the real functions £(z) and H(z) from the defin¬ 
ition, 

f -= 5 «<*>-in(z) , < 4 - n > 

o L-ikw 

and where $(K) is given by Eq. (2.7). Rather than work out the condition 
for validity of (4.10) in general, we will now specialize to the case 
of a focused beam observed at its focal point, i.e., we choose R = L. 

It is important to note that (4.10) ho^s only for locations on the cen¬ 
tral axis anyway. Let us introduce the focusing factor B f = kr Q /L, the 
ratio of the aperture radius to the radius of the diffraction-limited 
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(4.12) 


focal spot (8 f >> 1). It can be seen from (4.11) that 


£(z) = z/L 

H(z) = 6 f (L-z)/L , 


2 11/3 

Set <I>(K) = 32itYC K ' 
n 


in the focal plane, 
and substitute (4.12) into (4.10) to obtain 


where Y = 0.033n - 0.326 


^5x 2 ) = J c\ 2 J° dKK" 8 / 3 ^dze 


-8 f K 2 (L-z) 2 /kL 


jl-cos[K 2 z(L- z)/kL"| 


(4.13) 


Replace z by a new variable z' = z/L (we leave off the prime for convenience) 
and let x = R 2 (L-z)/kL. These two coordinate changes transform (4.13) 
into, 


= Y C n 2 k 7 / 6 L 11//6 t dz(l-z) 5 / 6 / i ”dxx‘ U / 6 £ 


-p (l-z)x 


4 


4 


[l-cos(zx)] 


= Yr(-5/o)C n 2 k' //6 L 11 / 6 

x f dz(1-z ) 5 / 6 {[8 f (l-z)] 5//6 - Re[p f (l-z) ± iz] 5 / 6 } 


This result is essentially equal to Eqs. (28) and (29) of Ishimaru 
[5], For large B f we utilize the asymptotic form for the braced factor- 
in (4.14), 


{—-> 


5_ 

72 


• 7 /Vd -z)- 7 '* 


} 


to obtain a beta function in integral form, which yields: 

•B f ' 7 / V k7/6i - u/6 - <4 - l5a> 

The numerical coefficient in Eq. (4.15) is 0.102. The focusing factor 
8^ being quite large in all cases of interest, it follows that ^^ 
is small compared to unity, and that is a sufficient condition for validity 
of the Rytov approximation. There is no practical need to work out an 
expression for <$x 2 > in the saturation regime because the log-amplitude 
variance appears to be small for all conceivable cases of interest. Let 
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(4.15b) 


us insert 8 - = kr /L into Eq. (4.15a) to obtain, 
f o ' 

< 5x 2 > = 0.102C 2 l\ * 7 / 3 , 


which is basically analogous to the geometrical-optics result foi ampli¬ 
tude fluctuations in a turbulent medium (see page 251 of ref. 2) with 
scale-size cutoff at l = T q . The result is not dependent upon wavelength. 

The conclusion that is ver y small for all practical situations 

(r - 0 . 5 m and L of the order of sever *1 kilometers) leads to an extremely 
important corolli-ry: amplitude fluctuatiom, are ^inportant in the foca l 
plane . It is therefore no coincidence that Eq. (3.3) can be obtained 
by a simple geometrical-optics ray-bending argument because refractive 
effects are dominant. To say it in yet another way, it appears sufficient 
to describe the distortion of the electric field on its way from the 


aperture to the focal pltae by simple geometrical-optics phase fluctuations. 

Let us now consider collimated beams, defined by allowing radius 
parameter R to become infinite in the definition of w in (4.6). In 


that case we find from the new definition and from (4.11) that 


# 2 4 

1 + zL/k r 

C(z) = -- 1 


n(z) = 


1 + L 2 /k 2 r 4 

(L-z)/kr 2 7 

- - - —~—7 = (L-z)/kr 2 

2/24 ' o 

1 + L /k r 

' n 


(4.16) 


When substituted into (4.10), n(z) clearly yields a negligible attenua¬ 
tion factor. Therefore (4.10) reduces to ti'e plane-wave form as we ex¬ 
pect, provided L « kr Q 2 (i.e., provided edge diffraction is negligible). 
Now consider diverging beams (R < 0) and particularly let 

R « kr 2 . It follows from (4.6) and from (4.11) that 
o 


i(z) = (z+R) (L+R)/(L+R ) 2 « (z+R)/(L+R) 


n(z) 


R . R(L-z) << ^ 
kr 2 (L+R ) 2 

O 


(4.17) 


Here again, n(z) yields a negligible attenuation in (4.10), and now 
( 4 . 10 ) reduces to the spherical-wave form for a point source radiatinr 
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2 


is not unduly restric- 


at the location (0,0,-R). The condition R « kr Q 
tive. 

In other words, the general beam-wave case is only of academic inter¬ 
est. Collimated and diverging beams can be treated, respectively, as 
plane and spherical waves with regard to amplitude scintillation. We 
have already performed a separate calculation for a focused beam which 
indicated that the .saturation effect is not observable in the focal plane. 
This entire treatment of beam waves appears to be as exhaustive as is 
required for amplitude effects. 

Some comments on phase are in order. In the model we have developed, 
there are diffractive and refractive effects. Diffractive effects are 
important for amplitude effects at all ranges; refractive effects start 
playing a role only in the saturation regime. With phase, refractive 
effects are important everywhere and it appears that diffractive ones 
can be ignored to a large extent. The most general expression we have 
for phase can be derived from t.q. (A4) of TRV, yielding 

6c(> - (k/8« 2 ) j** dzyd 2 K.6e'(K,z)e" lK ’ p(z) cosrK 2 (L-z)/2k] (4.18) 

o 

The last factor describes diffraction corrections which are important 
only from small eddies with size & < (I./k) ^ , but on the other hand the 
phase variance is determined chiefly by the large eddies with & >:> (L/k) 
Another way of showing the relative unimportance of the cosine factor 
in (4.18) is to determine the mean square of (4.18), ignoring the 
exp(-iK-p(z)) factor. After the usual universal approximations the co¬ 
sine factors reduce to 


cos 2 [K 2 (L-z)/2k] 


j £l + cos[K 2 (L-z)/k]^ 

1 - | jf 1-cos [K 2 (L-z )/k ] ^ 


(4.19) 


This shows that ^ (64>) 2 reduces to ^ (<5<}> g ) 2 ^ - K ? where <5<t> g is 

the phase (geometrical-optics) without the cosine factor. Our results 

indicate that K 6x/ « 1, nearly everywhere except in the initial part 

2 

of the saturation regime where 4 ^ can approach unity. Consequently 

/ / o ,\2 S S / r a . \^\ 













The significance of ray bending upon phase is easily illustrated 
from (A.18), ignoring the unimportant cosine factor, which can then be 
written, 

6<{) = kj ds6e(s), (A.20 

o 

where s is the ray-path parameter initiating at s = 0 and terminating 
at s (L). We compare 5(f) to 6tp g given by a straight path: 



(A.21) 


An estimate of the difference is given by 


6(Jr 6<t> =k J dz[ 6e(p(z),z)-6e(0,z)] 
8 o 

- f L - 

= kp (z) • J dzV T 6£:(0,z) 
o 

= kp(z)-[dp(z)/dz] 


„ I k r d P C£l] 

" 2 1 dz z=L 


(A.22) 


If we utilize the approximate relationship 

< p 2 (z)> - C n 2 < m l/3 z 3 , we find that 64>-64> g ~ . This is 

extremely interesting because the Rytov approximation has been shown 

to be valid when C 2 < ^kL 2 « 1. Thus the equivalence of and 
n m 

64> , the unimportance of ray bending, is once again demonstrated as a 
sufficient condition for the Rytov approximation. A corollary of this 
statement is that ray bending must be considered in the saturation regime 
where (A.20) can differ appreciably from (A.21). In short, Eq. (A.20) is 
a non-trivial extension of (A.21), the usual phase approximation, into 
the saturation regime. The question of the statistics of <54> cannot be 
fully answered at this time because there are not enough measurements in 
the saturation regime. Buser and Born [11] have observed discontinuous 
phase jumps that seem to point towards the crossing of diverse ray bundles 
in the image plane. More measurements are needed to confirm this phenomenon 
which has been questioned. 
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5. THE HERRMANN-BRADLEY NUMERICAL CALCULATION OF IRRADIANCE 


Equation (2.4) for plane waves can be written as 

I s + (2ik) _1 V 2 B - “• 6eB = 0 (5.1) 

0 z 1 ^ 


Our method of solving this equation for waves in turbulent air has been 
to convert it into an integral equation and to utilize small-angle approxi¬ 
mations. An alternative method has been to set B = exp'i' and to derive 
an equation for 4> from (5.1). Our saturation-regime results were obtained 
by taking ray bending into account in the equation for 4*, i.e., by in¬ 
tegrating along rays. Herrmann and Bradley [6] have attempted numerical 
integration of Eq. (5.1) using acceptable spatial-spectrum properties 
of 6e to select samples of an ensemble. They then obtain laser-beam pro¬ 
files of intensity that - ideally - should average to what can be cal¬ 
culated but have the advantage of giving as much detail as is desired 
(means and variances do not give much detail). A disadvantage of the 
method is that it costs time to produce a profile, and that it is perhaps 
difficult to produce enough to check the ensemble statistics. 

The main question is whether the numerical solution is valid or whether 
the approximations made in generating it are too severe. We address our¬ 
selves to this question here. 

Herrmann and Bradley convert (5.1) to an interaction representation 
by introducing the WKB phase increment 


54> g (p,z) dzf£ (£ Zl ) 


(5.2) 


and by transforming B to B = Bexp(-i6$). 
operator H, 


H 




Introduce a new differential 

(5.3) 


such that the Laplacian V 2 wcrks not only upon exp(i < $'t> ) but also upon 

A O 

whatever function of p follows H within a terminating parenthesis. The 
transformation B ^ "b and the definitions (5.2) and (5.3) then yield, when 


substituted into Eq. (5.1): 


3b/ 3Z - i(HB) = 0 


(5.4) 















Now, a formal but not very useful way of generating a full solution 
of (5.4) is to divide both terms by B and then integrate over z. This 
yields 


it 

B(z) = exp{-i y"dz 1 [B _1 (z 1 )HB(z 1 )] }. 


(5.5) 


The crucial approximation made by Herrmann and Bradley in solving (5.5) 

is to replace exp(x) by (l+x)/(l-x), where x is the exponent in (5.5). 

There are other approximations, but they are less critical. The ensuing 
* ... . & 

implicit approximate equation in B is then solved. If the approximation 
is to be a good one, we must demand that 


z 

|-i J dz^ [B ~ 1 (z )H$ , (z 1 ) ] | « 1 , 


(5.6) 


and the error is to third order in this parameter. Transform B back to 
■* ? 

B, and H back to to obtain 


£k/' i2 l B ' 1 <‘L ) V B < Z l>l « 1 . 


(5.7) 


Now replace 7 B(Zj) in (5.7) by the other terms of Eq. (5.1) to obtain: 


f Z 

2 ~J dz 1 <5e ( z 1 )^* 1 (z 1 )3B(z 1 )/3z 1 


(5.8) 


Replace B by B = exp(x + i<54>) where X is the log-amplitude and 6<j> 
the phase of the normalized field (both are real quantities), and utilize 
(5.2) to obtain 


ji [ 6 (fi - 6 c{) ] + x I « 1 

o 


(5.9) 


Certainly both real and imaginary parts of (5.9) must be small 
separately in order for the inequality to hold. It follows that X <<c 1 
is required. As X is a stochastic quantity, perhaps normal in distribu¬ 
tion (or at any rate not very differently distributed), it follows that 
its variance must be much less than unity. This restricts the 

Herrmann-Bradley procedure to situations where the Rytov approximation 














We have shown in section 4 that amplitude fluctuations are small 
and properly described by the Rytov approximation in the focal plane. 
Thus it follows that the Herrmann-Bradley procedure may yield accurate 
results for focused laser beams. On the other hand, the procedure may 
be seriously in error when utilized to describe collimated and diverging 

beams under conditions under which the parameter combination 

„ 2, 7/6 11/6 . 

C k ' L is not small, 
n 


20 






6. THE EFFECTS OF INTERMITTENCY 


Most of the fluctuations of laser-beam parameters such as irradiance 
and angle of arrival in turbulent air are directly due to fluctuations 
in the refractive index. A major portion of these are due to temperature 
fluctuations. The study of temperature fluctuations in turbulent air, 
e.g., as discussed by Tatarski [2 1 , yields the inputs to the analysis 
of these fluctuating phenomena. Parameters such as the refractive-index 

2 -i 

structure constant C and the scales of turbulence £ (or < ) and L 

o o m o 

derive from the study of temperature fluctuations in turbulent air. In 
recent years, workers at NOAA in Boulder, Colo. (R. S. Lawrence, G. Ochs, 
et al.) in particular, but not exclusively, have pointed out certain 
large-scale, low-frequency temperature fluctuations that have been ex¬ 
cluded in prior treatments. These give rise to a class of phenomena that 
do not appear to be sharply defined. Together with tne underlying cause 
these phenomena are referred to as "turbulence intermittency". 

Recently, Kerr [7] has discussed some of the effects of intermittency. 

He models the atmosphere into a number of irregular columns (slabs in 

2 

his specific treatment) with widely varying values of C^ , ind computes 
the effect upon the usual calculation of log-amplitude scintillation. 
Seemingly separate from this model is a discussion of finite-time averag¬ 
ing which can give rise to a data spread in apparently constant statis¬ 
tical parameters such as means and averages. 

It is not very clear, in our opinion, that the so-called intermittency 
of turbulence can be defined and treated as some phenomenon separate from 
what is otherwise just a random sampling of a statistical process with 
known probability distribution and parameters. We will try to approach 
a definition operationally as follows. 

Out of a host of statistical parameters describing the behavior of 
laser beams in turbulent air, certain variance parameters have been 
chosen for analysis because they are relatively fundamental and easy to 
measure. Two random variables, the variance of which we have analyzed, 
are the log-amplitude and the angle of arrival. A common mode of des¬ 
cription of both of these parameters as well as of others such as phase 
be giver in terms of integrals, 
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S(t) = Kr /dz. f d ? K f (K,z) 6 c (K,z,t) 
ATT 2 / 1 ^ 


( 6 . 1 ) 


Here, 6 e (K,z,t) is the transverse two dimensional Fourier transform of 
the deviation in dielectric permittivity 6 etr,t), and f(K,z) is a determi¬ 
nate function. For example, if F(t) is the log-amplitude, then f(K,z)= 

- (k/2) sin[ K (L-z)/2k] in the Rytov approxima:ion for plane waves. If 

£(t) is the angle of arrival vector, then f(k,z) is the vector -ik. In 

2 

order that £(t) describes the Rytov phase we set f(K,z)= (k/2) cos[K (L-z)/2k], 

■X 

For the purpose of this section it does not matter what form f(K,z) takes; 

Eq. (6.1) is sufficiently general for describing an operational definition 
of intermittency as it monifests itself in the optical phenomena. 

Now define a mean square C^. 2 = ^£ 2 ^ ; it is a variance as well as a mean 
square whence) =0. However, we cannot measure C^. 2 . 
is a quantity, 


What is really measured 


t-T/2 


( 6 . 2 ) 


If (6.1) is substituted into (6.2), then the time-dependent part to scruti¬ 
nize closely is 


t+T/2 „ _ 

dt' 6 e(K^,z^,t') 6 e(K 2 Z 2 ,t ? ) 

t-T/2 

Consider a random function g(t), with zero mean for convenience. 



(6.3) 


Let 


Mo-i / +I/2 d f g ( t ., . 

t-T/2 

and let <g> be the ensemble average of g. The temporal mean g^(t) will 
fluctuate around <gX and leas so as T increases. If the statistics are 
stationary, then <g> will be the limit of g^(t) as T -* 30 , and it will in¬ 
deed be independent of t. Hunt and Collins T 8 } have discussed this matter, 

2 

and it is rather easy to reconstruct their result for the variance a ^, of 
g T (t) around <g>. 
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o T 2 - (x/T)0 2 


(6.5) 


where a 2 is the ensemble variance of g(t), and T is the integral time of 
the normalized auto-correlation of g(t). Eq. (6.5) is correct only for 
T << T. In the opposite limit, will vary randomly as g (t) does be- 

■t ^ 

cause the frequency 2nT lies beyond the frequencies of interest in the 
power spectrum of g(t). When T <<T , we may replace the time average by 
the ensemble average to good approximation. 

Consider the time T Q to be the longer of L q /U t or L q / 6U where U T is 
the wind velocity component across the laser beam, and <SU is the (more or 
less) isotropic root-mean-square wind velocity. Let averaging time T be 
much greater than T . Our cardinal assumption for (6.3) is 


1 

T 


^+T/2 

J 

t-T/2 


dt , 6^(K 1> z 1 ,t')5e(K 2 ,z 2 t’) 


( 6 . 6 ) 


-i r 

T t-T/2 




'' t ' C n (Z l’ t,)$ 2 (K,Az)6(K l +K 2 ) 


where <t> 2 (K,Az) is the partial Fourier transform of the autocovariance of 

6e(r) denoted as F(K x> K y ,Az) by Tatarski [2], et al., and 5 ( K 1 +K 2 ) is a 

two-dimensional Dirac delta function. The assumption implies als> that 

the time variation in C 2 (z,t ) is slow compared to the time T q defined 

n 

above. We now insert (6.6) into Eq. (6.2) and utilize the universal 
approximations that center around the rapid decay of $ 2 (K,Az) with in¬ 


creasing Az to obtain 


C 2 (. t , T )=-^2 J iz J d2R f2 ( K ’ z) $(K) 


t+T/2 , 


1 • 2 ’ 

X f dt C„(Z,t ) 


t-T/2 


(6.7) 
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where $ (K ) is given by (2.6) and (2.7). The variance C^(t,T) is a 
fluctuating quantity to the extent that the last integral factor is. 
C 2 (z,t ) is time dependent, or - weaker statement- if the integral 


“ 9 • 9 -* 

Jdz C n Z (z,t ) f Z (K, : 


is a time varying quantity, then so is C^(t,T). Considering it established 
that T >> T , we can see two possible causes of time variation: 

(i) hon-stationary statistical behavior of C n 2 (z,t ) because T is so 
long that the heat input into the atmosphere changes (diurnal 
patterns). 

(ii) Large-scale spatial variation in C n (z,t ), associated with lengths 

L such that time L /U is not much smaller than T (obviously L c » L q ). 

c c # 

This second effect appears to be most basic to the concept of intermittency. 

In order to proceed more specifically, we assume a large-scale frozen- 

flow model. That is to say, we assume that large-scale variations in 

C 2 (z,t) are blown across the line of sight by a wind with velocity vector 

n /v 

U (unit vector U) not necessarily perpendicular to the line of sight. The 
assumption implies at point r= (0,y,z): 


9 ^ * 94-*' 

C n 2 (r,t )= C n (r- U(t -t), t) 


( 6 . 8 ) 


We will utilize this in (6.7). However, first we will assume a Kolmogorov 
spectrum derived from Eq. (2.7) for *(K) in (6.7) and formally perform the 
d 2 K integration. Thus (6.7) becomes 


, ^ 1 t+T/2 2 
I(t,T)= / dz f(z) - f dt' c n ) 

• / o T J 

t-T/2 


(6.9) 


The function f(z) has the following form: For log-amplitude, f(z) 

0.563 k 7/6 (L-z) 5/6 . For phase it is proportional to k 2 L Q , and for 
the focal-area parameter r it is 3.9k (L-z) . Similar and related 

Lb in 
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forms for f(z) result, depending on which parameter £(t) is studied in 
Eq. (6.1). Now insert (6.8) into (6.9) and transform time variable t to 
spatial coordinate s * Ut (s measures distance in the wind direction): 

L UT/2 


C r 2 (t,T) = J dzf(z) dsC^ 2 (r-sU,t) (6.10) 


n 


The mean of (6.10) yields, under the assumption (for simplicity) that 
2 

the mean of C does not depend on r. 


L, 

<Cf 2 > - / d*f(z)<C n 2 > 


( 6 . 11 ) 


In order to compute the variance of C r (t,T), we will introduce a spec- 

•* ^ 2 
trum $ (K) of the large-scale fluctuations of C (r) via the transform 
c n 


<V <'>„ «,» - W > 

■ —J’ 1 

8tt j 


( 6 . 12 ) 


yd (K)exp (-iK» ( r 1 _ r 2 ) ] . 


Using (6.12), we can compute the variance of (t,T) from the average 
of the square of (6.10) minus the square of (6.11). The two resulting 
ds integrations are easily performed to obtain, 


<V> - <c £ 2 ; 2 


-\/d 2 K* <K) 
4tt c 


(6.13) 


x I Si g?lr/2 /2) f dzf(z)exp(iK,r)| 
o 

The factor |sin(K*UT/2)/(K*UT/2) | 2 acts as a low-pass filter particularly 

in directions normal to ”z, upon the spectrum $ (K) of fluctuations in 

2 c 
C n (z). It suppresses contributions from scale sizes 1 « U^T/2 in 

analogy to the time-dependent case in (6.5). 

Intermittency thus causes data spread in the variance parameter for 

relatively short averaging times T, and the variance of this data spread 

2 

for a seemingly steady parameter (which would be constant at very 
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large T under statistically homogeneous conditions) is given by Eq. 

(6.13). Consider as an example the log-amplitude £ = InA in which case 

C „ 2 i s the log-amplitude variance. Then, Eq. (6.13) expresses the 

variance of the data spread of a finite-time log-amplitude variance in 

terms of the wind velocity U and the large-scale fluctuations of 

C 2 (z) that constitute 4> (K). We not? that the data-spread variance de- 
n c 

creases with increasing wind velocity when U^,T/2 is sufficiently large 

with respect to the important scale size.?. This trend is clearly visible 

in the log-amplitude variance measurements by Kleen and Ochs [9]. 

Intermittency, then, appears to be a large-scale variation in 

C 2 (z) as far as the optical effects are concerned. It affects the sta- 
n 

tistics of optical quantities such as phase and amplitude when these sta¬ 
tistics are determined over time periods not very long compared to the 
time in which a typical large-scale variation crosses the beam. Its char¬ 
acteristics are determined by the (unknown) spectrum ^(K) which in all 
likelihood is highly anisotropic and strongly characterized by local 
terrain and weather conditions. 

The large-scale variations in refractive index probably have no di¬ 
rect optical consequences because both refractive and diffractive bending 

re governed by small-scale variations. An intermittent variation of 

2 2 

.tatistical parameters such as and is brought about. Perhaps 

this variation should be regarded in the same way that a short-term (1-5 
seconds) measurement is an instantaneous realization of which one can 
only give the statistics. Here, that would imply knowledge of (6.13) 
which in turn requires knowledge of ^(K) which is probably not a uni¬ 
versal function but a strong one of terrain and weather conditions xocally. 
That may render attempts to determine $ c (K) in general ineffective a 
priori. In this sense it remains unclear whether intermittency phenomena 
can be distinguished from random sampling of a statistical process with 
known probability distribution. 










7. ANGLE-OF-ARRIVAL SPECTRUM AND MICROSCALE 


The power spectrum of the angle of arrival of a ray has been cal¬ 
culated in TRIV, and an extension of this calculation for the angle-of- 
arrival difference of two rays is given in TRVI. An important aspect 
of this power spectrum (we consider it for one ray because the inter¬ 
ferometer situation does not change matters significantly) is its high- 
frequency behavior. Let (b = U t /L q and ^ (we assume frozen flow 

for the time being; the alternative situation has been discussed before). 
The power spectrum Wq(gj) can then be written for u) >> as 


W 0 (o>) = 4yC 2 LL _1/3 W fi (o>,5/6), 
t) n o o 

„ / c/n rq/3 ) /n /W_n~ 2/ 

W Q (ai,5/6) - r(5/6) ^ <^> 


for w < 


(7.1) 


n -2/3 -1 -(w/JLr 

<ir> • 

T T T 


-2/3 

That is to say, there is a very sharp transition from an <a dependence 

2 2 

for u) < to a steep gaussian decay ~ exp(-w /fyp ) when w exceeds 

The question arises whether this steep transition can be used to obtain 

a microscale measurement by determining under steady wind conditions 

(so that < may be retrieved from ft = < U T ). 
v m T m T 

It appears to us that non-zero aperture limitations prevent an ob¬ 
server from finding this transition, even in the idealized noise-free 
case. In actual fact, the beam is not infinitesimally thin; it has a 
nor-zero thickness d (e.g., d is a gaussian halfwidth parameter of a 
thin collimated beam). Temperature fluctuations at scale sizes i larger 
than d do indeed displace the entire beam (beam wander) as if it were a 
ray, but smaller ones with l < d decompose the beam into parts (beam 
spreading). It is quite true that the sum of both effects would be 
observed in long-term statistics to be governed by a variance 

Q 2 > ~ c 2 L£ -1/ ^ 3 where i is the microscale, but the measurement is 
no o 

of the displacement of the center of irradiance, or of some related 

i 2 v 2 -1/3 

quantity which is governed by a variance \9 / C n Ld * Tatarski [2] 







has computed a power spectrum for beam wander of a beam of thickness d, 

and he finds a spectrum that resembles (7.1) for u) < U/d. However at 
u) ~ U/d the minus two-thirds behavior is terminated and a transition to 
a minus eight-thirds power of frequency occurs. The transition point 
is now characteristic of d rather than of £. A second transition at 
u) ~ ft is no longer sharp and presumably hidden in noise. 

The spectrum of (7.1) would indeed yield microscale information if 
aperture d were less than Unfortunately can be very small (milli¬ 

meter or less) indeed and the intensity of such narrow beams are probably 
too weak for use in this way. Therefore the limitations in extracting 
a microscale measurement from an angle-or-arrival spectrum appear too 
stringent at present. 
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8. MEAN FLUCTUATION RATES 

The power spectrum of a fluctuating quantity contains much information 
about the time scales or the fluctuation rates of that quantity. For 
example, the power spectrum of a fluctuating amplitude of an electric field 
tells the energy of the field as a function of frequency. The integral 
of the power spectrum over all frequencies is simply the total energy (on 
the average) of the fluctuating field. A range of frequencies (fluctuation 
rates) ranging from U/L q to U/£ q roughly (where U is a typical velocity 
driving spatial irregularities across the laser beam) contributes to this 
total energy. One might wonder whether there is a dominant frequency. 
Perhaps inspection of the power spectrum might show up such a dominant 
frequency. However, it requires quite a bit of data processing to produce 
a power spectrum, and then it might be somewhat difficult to read a 
dominant frequency out of it. Furthermore, it is convenient to lave a 
simpler more direct measure of the rate at which a random quantity, say 
the irradiance, fluctuates. A simple measure is the average number of 
zero crossings. The observer takes a certain time record of the irradi¬ 
ance, determines a mean irradiance level during that time, and then counts 
the number of times that the actual time record of irradiance crosses the 
mean level. It is clear, intuitively, that this count is a measure of the 
graininess, the irregularity of the fluctuating quantity. 

Beckmann [10] gives the mean frequency of zero crossings as o f /7r where 
the variance parameter of the random quantity f is given bv 


a. 


= iL r<f(t)f(t+r)> 

' - 3t 2 L<f(t)f(t)> j T . 0 


R f "(0) 

R f (0) 


(3.1) 


We will compute this quantity for angle of arrival, amplitude, and phase. 

Angle of Arrival - The autocorrelation cf the angle of arrival of a 
ray is given by Eq. (31) of TRIV which is written equivalently as. 


R e (T) = 2yC, > 


2. -1/3 


/ 


. 2 t 2 
, - ,, -x/K L 
, ... .-11/6 m o 

dx x(l-hx) e 


( 8 . 2 ) 


X J o ( V ^)e- ,!<i “ T/2) ' 
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where y = 0.033TT 2 ~ 0.326, ^ - U t /L q , = 4 AU/L q /3. 

This autocorrelation yields the variance for T = 0, namely 
R (0) - 3.9ir(7/6)C 2 k 1/3 L. The second derivative of the autocorrela¬ 
tion at x = 0, namely R"(0) is easily computed from the second line in 
Eq. (8.2) which yields 




[j < V /x>e-* (W2 > 2 ] 

L ° T T=0 


= - § ( 0 ) W) = - f(v"/L o i ) 


(8.3) 


where v 2 = U 2 + 16AU 2 /3 is a variance of velocity that combines the 
frozen-flow and the random-flow velocities into a single parameter. Un¬ 
less U « 0, the difference between v and U T is negligible. Using (8.3) 
in (8.2) we apply the definition (8.1) for angle of arrival (f = 9) to ob 

tain 



r 
21 . 2 J 


2 —11/6 ' x/,C m 

dxx (1+x) e 


V//' 


- 11/6 -* lK m \ 2 
dxx(1+x) ' e 

(8.4) 


In numerator as well as in denominator, it is permissible to replace 
1+x by x because » 1. It is then easily established that 


°e 2 - f ( v/ V 2 


(8.5) 


Thus a non-surprising result has been found in yet one other way: the 
angle-of-arrival fluctuations are governed by the smallest turbulence- 
scale sizes. In practice (see section 7) the graininess will probably 
not be determined by i Q but bv the diameter d of a pencil beam by which 
a ray is simulated. Thus, the fluctuation rate vM q or v/d is a measure 
of the most important rate of angular beam deviations. 

Loe-amplitude - The log-amplitude variance in the Rytov approximation 

is given by 


<«X 2 > - YC n V / 6 I . 11/6 /°to' 11 / 6 <l-x^sin*) 


( 8 . 6 ) 
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in which expression the variable x was obtained from the wavenumber vari- 

2 

able K through the transformation x = K L/k. As usual, 

Y = 0.033^ 2 - 0.326. We generalize <^5x 2 ^ = < <$X(t)6x(t) > to 

R a (t) = <5x(t+T)6x(t)> by the same procedure as in the case of angle- 

of-arrival, i.e., we replace <I>(K) by 


4>(K)J Q (KU T T)exp(- | KAU 2 t 2 ) , (8.7) 

(see Eq. (30) of TRIV). After transforming K to x, and utilizing the fre¬ 
quency variables oj and Aoi, we insert the extra factors in (8.7) into 
(8.6) to obtain 


•00 ™l/y /T 

■ m-ycW 1 ' 1 / " 

A n J 

o 

T , -y(A(jJl/2) 2 

x J o (a) T i/y)e 

where y 5 x(kL Q 2 /L). The calculation of a 2 = -R^"(0)/(0) 
repeated by utilizing the intermediate step (8.3) to obtain, 


( 8 . 8 ) 


is easily 



kv 

2L 


f 


r/t J “kx/K 

, ”5/6,, sinx N m 

dxx (1- ——)e 


7 f' 


, -kx/< L 

, -11/6., sinx. m 

dxx (1--—)e 


(8.9) 


The integral in the denominator can be approximated well by ignoring the 

exponential factor. It then yields the factor n/[2r(17/6)sin(7TT/12) ] in 

the denominator. The integral in the numerator can be approximated by 

2 1 / 6 

ignoring the sinx/x term, and it yields a factor r(l/6)(< m L/k) . The 
result is 

/< 2 L 

2 _ r(l/6)r(17/6)sin(7TT/12) m 

a A " TT \ k 

- 2.95(k 2 /k) 1/6 (kv 2 /L) 

IQ 

2 1 /6 

Bearing in mind that the factor (< L/k) is of order unity in 
° m 

optical propagation problems of interest, we note that the rate of zero 

1/2 

crossings is given by the frequency v(k/L) , i.e., by those irregularities 

i /o 

with scale size (L/k) . That is in good agreement of course with the 


\l/6 
/ < 


(kv 2 /L) 


( 8 . 10 ) 
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interpretation of the spatial covariance of log-amplitude which is also 

1/2 

governed by a length - (L/k) 

Phase - Finally, the calculation is easily repeated for the pnase in¬ 
crement which we write 


6$ 




dz6e(z,t) , 


( 8 . 11 ) 


where a somewhat unimportant filter factor has been ignored. Define 
r + (t) * <<54>(t)64)(t+T)^ • We obtain it from R^(0) by means of the sub- 


4> 

stitution (8.7) and get 


-K% 


R^(T) * AyC 


2.. -2. -11/6 vn id.. 2 


dKK(K +L q ) e 


2yC n 2 k 2 LL o 5/3 ^° dx(l+x) -11/6 J o (a) T iv / x)e 


J (KU_T)exp[-A(KAU t) /3] 
0 T 


-x[ (Acot/2) 2 + 1/k 2 L 2 ] 

mo 


( 8 . 12 ) 


2 

Note that this yields the well-known phase variance R, (0) = 64> = 


(12/5)YC n 2 k 2 LL Q 5/3 . Utilizing (8.3) we obtain 


R,"(0) = yC 2 k 2 LL 

4) n 




/ 2 t 2 
,. -x/< L 

. -v —11/6 m o 

dxx(l+x) e 


(8.13a) 


which can be approximated well by replacing fl+x) >y x to yield, 


1/3 2 2 -1/3 2 

V«» -Yr(l/6)0C B L o ) /3 C n kLL 0 V 


(8.13b) 


Thus we find for the zero-crossing rate parameter, 

°4, 2 ■T2 r(l/6)(,: m 1/3v2/L o 5/3) • (8 ' U> 

which result indicates that phase fluctuations appear to be governed by 

a rate ~ v/(L ^ 6 ), which is not very different from v/L because 

1/ fi 0 0 0 
(l /L ) is of order unity, 
o o 

The calculations of zero crossings thus seem to indicate that angle 
of arrival, amplitude, and phase each have their jwn typical fluctuation 
rate f namely, 
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angle of arrival : f n ~ v/i. 

° ° 1/2 

log-amplitude : f ~ v/(L/k) 

phase : f ~ v/L 5/6 i l 1/6 - v/L 

4> O O 0 

The full power spectrum indicates a spread of rates between v/L and 

o 

V /*o (there is an artificial smoothing of the spectrum between to = 0 and 
w ~ v /L o th’t has no physical meaning) but these typical rates are dominant 
in the sense that the average number of surges and fades of signal records 
are given by them. 
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